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RESOLVED SHOCK STRUCTURE OF THE BALMER-DOMINATED FILAMENTS IN TYCHO'S SUPERNOVA 

REMNANT : COSMIC RAY PRECURSOR? 
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ABSTRACT 

We report on the results from Ha imaging observations of the eastern Hmb of Tycho's supernova 
remnant (SN1572) using the Wide Field Planetary Camera-2 on the Hubble Space Telescope. We 
resolve the detailed structure of the fast, collisionless shock wave into a delicate structure of nearly 
edge-on filaments. We find a gradual increase of Ha intensity just ahead of the shock front, which 
we interpret as emission from the thin {r^ 1") shock precursor. We find that a significant amount 
of the Ha emission comes from the precursor and that this could affect the amount of temperature 
equilibration derived from the observed flux ratio of the broad and narrow Ha components. The 
observed Ha emission profiles are fit using simple precursor models, and we discuss the relevant 
parameters. We suggest that the precursor is likely due to cosmic rays and discuss the efficiency of 
cosmic ray acceleration at this position. 
Subject headings: ISM: individual objects (G120.1+1.4) — ISM: supernova remnants — shock waves 



1. INTRODUCTION 

The shock transition in fast astrophysical shocks is in- 
trinsically a "collisionless" process, and energy is dis- 
sipated via plasma turbulence and/or electromagnetic 
fields. An important consequence of the collisionless 
nature of the shocks is cosmic ray acceleration (e.g., 
iBlandford fc Ei chlcr 1987). While there is increasing ev- 
idence of cosmic ray acceleration in supernova remnants, 
the details of the process are still not well understood, 
and the question of whether supernova remnants (SNRs) 
are the prim ary ac celeration sites of Galactic cosmic rays 
is still open ()Buttl l2009). 

Cosmic ray acceleration models require a precursor 
in which accelerated particles can scatter back to the 
postshock region for further acceleration. Observations 
of the cosmic ray precursor can constrain the two key 
parameters of acceleration models; the diffusion coeffi- 
cient and the injection effic iency (jBlandford fc Eichleii 
Il987t Boularc sfc Coxl 119880 . The Balmer-dominated 
filaments that are produced when fast SNR shocks 
propagate into partially neutral gas are potential sites 
where such a cosmic ray precursors can be observed. 
Most of the Balmer emission comes from a very nar- 
row zone behind the shock, where the hydrogen atoms 
swept up by the shock are excited befor e they are ion - 
ized ( Chevalierfc Ravmondl il978: Chevalier et al.lll"980l ). 
However, in the presence of the precursor, additional 
Balmer emission is expected from the precursor region 
where the preshock gas is compressed and heated. Using 
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long-slit Ha spectroscopy along the shock normal of a 
Balmer filament in Tycho's SNR. lLee et~all (|2007L Lee07 
hereafter) found that there is an increase of the Ha nar- 
row component intensity in a small region (^^ 0.4") ahead 
of the shock front, which they proposed as potential emis- 
sion from the precursor. However, the angular resolution 
of the observation by Lee07 is ~ 0.5" and their results 
needed to be verified with high resolution observations. 
In this Letter, we report Ha imaging observations of 
Balmer-dominated filaments in Tycho's SNR using the 
Wide Field and Planetary Camera 2 (WFPC2) on the 
Hubble Space Telescope (HST), which resolves the de- 
tailed structure of the shock. § [2] presents the observa- 
tions and reports the detection of the precursor. The 
precursor is modeled in § [3] and its characteristics are 
discussed in § |4l Finally in § [5l we discuss the efficiency 
of cosmic ray acceleration in this region. 

2. OBSERVATIONS AND RESULTS 
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Fig. 1. — Hubble Space Telescope image of the Balmer- 
dominated filaments in the northeastern part of Tycho's SNR. The 
WFPC2 detector with F656N filter was used. The inset shows 
an Ho image toward Tycho's SNR from a ground-based telescope 
(KPNO). The area observed by HST is marked. 
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We observed the eastern lirab of Tycho's SNR us- 
ing WFPC2 on HST. The observations were conducted 
on Mar 23, 2008, with the camera arranged to cover 
knot g, one of the brightest Ha filaments in the rem- 
nant (^2000, ^2000=00^ 25"^ 52% +64° 09' 21"). A to- 
tal of 10 exposures, each with 2600 ^ 2700 sec, were 
obtained using the F656N Ha filter, which has a cen- 
tral wavelength of 6564A. The filter has a bandwidth of 
54 A (~ 2500 km s^^), which is comparable to the ve- 
locity width of the Ha broad component in this region 
(jGhavamian et aLll2001l ). and transmits about half of the 
broad component flux. Images are combined and "driz- 
zled" onto a 0.06" pixel~^ scale using the IRAF task 
multidrizzle (jFruchter fc Ho ok 2002), which also detects 
and removes cosmic-rays. The PSF of the final drizzled 
image has a FWHM of ~ 0.18". 

The full field observed with HST is shown in Fig. [TJ 
The brightest filaments, comprising the knot g region and 
associated filaments, are seen on the WF3 chip, while 
fainter filaments belonging to the northeastern limb of 
the SNR are visible on the WF4 chip. The superb an- 
gular resolution of HST reveals details of the filaments 
not available from ground-based telescopes. This is more 
clearly seen in Fig. [2j where we compare the close-up 
view of the knot g region to the image observed with 
a ground-based telescope (CCD image obtained at the 
KPNO 2.1m telescope on 2007 October 4) at the same 
scale. The HST image reveals a faint extension of the 
emission toward upstream (to the east) along most of the 
bright part of the filaments (between (5=64° 08' 40" and 
^=64° 09' 20"). This extension is more clearly demon- 
strated by the Ha brightness profiles from the cuts along 
the shock normals, as seen in Fig. [3] The profiles show 
a bright emission peak of thickness < 0.5", which is the 
emission from neutral hydrogen excited in the postshock 
area. The small bumps around offsets —1.8" and 0.7" 
are likely due to the projection of fainter tangencies of 
the rippled shock front to the line of sight. The Ha 
emission slowly falls off not only toward the downstream 
but also toward the upstream direction. In the down- 
stream, the Ha emission is emitted within a very narrow 
region behind the shock as neutral hydrogen is rapidly 
ionized. Thus, we consider that the downstream emission 
is likely a projection of the curved shock fronts. Some of 
the upstream emission could be in principle attributed 
to a similar projection effect. However, Fig. [2] shows 
that the shock front in this region does not show any 
significant curvature, and it is difficult to devise a shock 
geometry that explains both upstream and downstream 
emission. Fig. |4] shows emission profiles along different 
cuts, and the upstream emission components are simi- 
larly seen. The upstream emission is also seen faintly in 
the profile of cut 01 where the shock surface is appar- 
ently convex (see Fig. [2]) , and a projection effect is not 
likely to explain the upstream emission. Therefore, we 
propose that the faint emission from the upstream region 
represents the emission mostly from the neutral hydro- 
gen atoms excited in a shock precursor, and the effect of 
shock geometry makes only a minor contribution. 

3. MODELING THE PRECURSOR EMISSION 

To estimate the physical properties of the precursor, we 
model the observed spatial Ha emission profiles (Fig.|4]). 
We first try a simple toy precursor model. We assume 
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Fig. 2. — (a) The Ha image of knot g, one of the brightest 
Balnier-doniinated filaments in Tycho, taken by KPNO 2.1m tele- 
scope, (b) The HST image of the same region. The black rectangles 
are regions from which the brightness profiles are extracted. 

a shock front located at a; = xq with increasing x to- 
ward downstream. The Ha emission from the postshock 
region is assumed to be confined to a region of thick- 
ness of w with a uniform emissivity. In the precursor, 
the Ha emissivity peaks at the shock front and decreases 
exponentially away from the shock front. The profile is 
written as 




xa)/L 



li X — xq > w 
\i Q < X — xq < w 
\i X — xq < 



where F is the total Ha fiux from the postshock region, 
L is the precursor length scale, and / is a flux ratio of 
the precursor to the postshock region. The variation in 
the downstream region could be modeled with a curved 
geometry of the shock. However, we find that the ob- 
served variation cannot be adequately fit by shocks of 
a simple geometry. Instead of introducing arbitrary fit 
parameters to describe the structure along the line of 
sight, we simply assume plane- parallel shocks and ignore 
the data where the model deviates from the observation. 
As long as L is sufficiently smaller than the local cur- 
vature radii of the shocks, the plane-parallel assumption 
will not significantly affect estimated precursor parame- 
ters. As evident from Fig. [51 the profiles require multiple 
shock components projected along the line of sight. To 
minimize the number of free parameters, we assume that 
shocks have the same profile shape, i.e., parameters /, L 
and w are tied among multiple shocks for a given pro- 
file cut from Fig. [2l The model profiles are Gaussian 
smoothed to account for the instrumental profile. 

The fits are shown in Fig. 31 and the results are sum- 
marized in Table 1. We find that the Ha flux from 
the precursor region is comparable to that of the post- 
shock area (f ~ 1), with L around 1", corresponding to 
3 X 1 0^^ cm at the assumed distance of 2.1 kpc to Ty- 
cho (iKamper fc van den Berghl [19781 JGhavamian et all 
l2001f ). The thicknesses of the postshock emitting area 
are around ~ 0.3". A relatively large L of ~ 2" is found 
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TABLE 1 

Fit parameters from the toy precursor model 

1 Tn ^] xVd.o.f. 



cut 01 


1 4+0'3 
-■-■^-0.2 


9 c: + 0.7 


cut 02 


n 79 + 0.06 
'J- '^-0.06 


1 n+0.1 
J^-^-0.1 


cut 03 


PI c:7+0.04 


0.70t°;°« 


cut 04 


„ 0^+0.04 

U.»0_g 03 


0.87«:«^ 



0.42 
0.31 
0.35 
0.39 



+0.06 

0.05 
+0.02 

0.02 
+0.01 

0.01 
+0.02 

0.02 



206/170 
170/182 
258/182 
220/182 



in cut 01, but this is likely to be overestimated due to the 
much stronger local curvature of the shock in this region. 
We also note that, for profiles in cuts 02 and 03, where 
the structure of the overlapping shocks is quite complex, 
the fitted precursor parameters are somewhat sensitive 
to the assumed baseline and also the number of shock 
components. These can cause '--^ 20 — 30% systematic 
uncertainties on the fitted parameters. 

We now consider a more realistic precursor model that 
assumes an exponential temperature profile (similar to 
the above toy model) in the precursor. The model calcu- 
lates the ionization of hydrogen atoms throughout the 
shock and the emissivity of Ha and Ly/3 lines. The 
radiative transfer of the Ly/3 line is computed using 
the Monte Carlo technique to account for the Ha en- 
hanceme nt by Ly/3-trapping. The model has been uti- 
lized by iWagner et al.l ()2009[ ) to interpre t the result of 
Lee07. We adopt the parameters used in IWagner et al.l 
(|2009f l: shock velocity of 2000 kms~^, preshock density 
of 1 cm"*^, and the preshock neutral fraction of 0.85. 
More detail s of the model and the input parameters can 
be found in IWagner et a l.l (!2009). The extensive discus- 
sion of the detailed modeling and associated uncertain- 
ties is beyond the scope of this Letter, and here we simply 
present a brief summary of the results. For cuts 02, 03, 
and 04, we estimate peak temperatures in the precursor 
Tpoak = 80, 000 - 100, 000 K, and the precursor length 
scale L = 5 ^ 7 X 10^^ cm. As in the toy model, a larger 
length scale is required for cut 01. The estimated length 
scales of the precursor are generally larger than those es- 
timated from the simple toy model. This is because the 
Ha emissivity is sensitive to the temperature, i.e., the 
emissivity profile increases more rapidly than the temper- 
ature profile, thus effectively reducing the length scale. 
We note that the precursor para r neters are in agreement 
with the results of lWagner et alj (|2009[ ). while the length 
scale is slightly larger. 

4. ORIGIN OF THE PRECURSOR 

The existence of a thin precursor has been suggested 
from previous observations. The spectral profiles of the 
Ha narrow component traces the velocity distribution of 
the gas entering the shocks, and those observed in SNRs 
have widths of 40 - 6 kms"^ (e.g., iSmith et al.l 119941 : 
iSollerman et al.ll2003() . This is too large for the temper- 
ature of the ambient gas, as all the hydrogen would have 
been ionized at the implied temperature and no Balmer 
filament should exist. Instead, the observed line width is 
suggested to represent the gas heated in the precursor, 
which is thin enough for the preshock neutrals not to be 
completely ionized. Also, the observed flux ratio of the 
Ha broad component and the narrow component was 
sometimes found t o be sm aller than what models pre- 
dict (e.g, .Ghava mian et al.ir2003; .Rak owski et al.i ,2009,) . 



/I 


precursor 

1 

1 


r- postshock 


^ Shock Direction 









-8 -6-4-2 2 4 6 

Distance from the peak f"] 

Fig. 3. — H« brightness profile across the shock normal extracted 
from cut 04. The bright narrow peak is believed to correspond to 
the emission from the narrow region of the immediate postshock 
area. The emission extends up to 4" toward upstream, which is in- 
terpreted as emission from the shock precursor. The small bumps 
in the profile (marked with small vertical bars) are likely the pro- 
jection of other shock fronts. 





Fig. 4. — Ha brightness profiles across the shock normals in dif- 
ferent regions marked in Fig. [2] The profiles are fitted with multiple 
components of toy precursor models. The gray area represents the 
accumulated emission from all models, and the solid lines represent 
emission from projected individual shocks. 

and the excessive narrow component emission was at- 
tributed to the contribution of emission from theprecur- 
sor (Ghavamian et al. 2001: Rakowski et aLll2009[ ). 

The characteristics of the precursor revealed by our 
HST observations are consistent with results from previ- 
ous observations. The peak temperature in the precursor 
may be relatively higher than the temperature implied 
by the line width of the Ha narrow component. How- 
ever, the temperature we modeled is the electron tem- 
perature, which might not be in equilibrium with the 
neutral hydrogen atoms. Also, the velocity profiles of 
the Ha emission could deviate from a Gaussian profile 
(Ghavamian et al. 2001; Lee et al. 2007; Ravmond et alJ 
[2010f l. and a simple line width may not be an adequate 
temperature indicator. Our observations show that the 
emission from the precursor is a significant contributor 
to the Ha narrow component. We estimate that the 
precursor emission may contribute up to 30-40% of the 
narrow component for a slit width of 1". The narrow 
component emission from the precursor will affect esti- 
mates of electron-ion equilibratio n based on the observed 
broad-to-narrow inten s ity ratio (jGhavamian et al.ll2001l : 
iHeng fcMcCravl [20071 : iRakowski et al.ll2009[ ). For those 
SNRs where the observed Ha broad-to-narrow flux ratio 
was smaller than the model predictions, accounting for 
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the precursor could bring them into agreement. 

The hkely candidates for producing this precursor 
are a cosmic ray prec ursor or a fast neutral precursor 
(ISmithet all [1994: Hes ter et al.l[T99l . While our HST 
results provide detailed structure of the precursor, mod- 
els of these precursors are not available for a quantitative 
comparison. However, the cosmic ray precursor scenario 
has been preferred over the fast neutral precursor for 
various reasons (see Lee07 and references therein). Also, 
the growing evidence of cosmic ray acceleration in Ty- 
cho's SNR supports the e xistence of cosmic ray precursor 
C e g.. I Warren e t al. 2005). 

I Wagner et al.i (,2009. ) computed a series of time de- 
pendent numerical simulations of cosmic ray modi- 
fied shocks, tuning the model parameters to reproduce 
the precursor chara cteristics of Lee07. The results of 
IWagner et al.l (2009) can be a plausible approximation 
as the estimated precursor properties are not much dif- 
ferent. They found that, assuming a distance of 2.1 kpc 
to Tycho's SNR, the CR diffusion coefficient, k, the in- 
jection parameter, e, and the timescale for the energy 
transfer, t, of k = 2 x 10^** cm^ s~^, e = 4.2 x 10"'^, and 
r = 426 yr is required to describe the observations. The 
length scales of the precursor estimated from our new 
observation are slightly larger than that of Lee07. The 
larger precursor length scale requires a larger cosmic ray 
diffusion coefficient. This increases the acceleration time 
scale, so the cosmic ray injection parameter may need 
to be increased to compensate for the slower accelera- 
tion. While some fine tun ing of the par ameters may be 
required, we believe that IWagner et all 's findings, e.g., 
that the cosmic ray acceleration is not very efficient in 
this shock and about 10% of the shock energy has con- 
verted to cosmic rays, are still valid. 



5. COSMIC RAY ION ACCELERATION IN TYCHO 

The existence of the cosmic ray precursor does not nec- 
essarily imply efficient cosmic ray acceleration. For the 
Balmer-dominated filaments to be observable, some neu- 
tral hydrogen needs to survive ionization in the precur- 
sor. As efficient cosmic ray acceleration tends to make 
the precursor wider and hotter, the Balmer-dominated 
filaments may not trace t he shocks having effi cient cos- 
mic ray acceleration (cf . iHelder et al.l l2009f ) . For the 
shock studied in this paper, a more direct suggestion of 
inefficient acceleration can be inferred from the results 
of LeeOT. From the difference in radial velocity of the 
preshock gas and the Ha narrow component emitted in 
the postshock region, they estimated the amount of gas 
deceleration in the precursor. The preshock gas is decel- 
erated in the precursor due to the gradient of the cosmic 



ray pressure by about a few hundred km s^ based on the 
radial velocity measurements. The value is not sensitive 
to the assumed distance, and is significantly smaller than 
the measured line width of the broad component. In this 
shock, the thermal pressure of the ordinary gas still dom- 
inates over the cosmic ray pressure and the acceleration 
is not likely to be very efficient. 

Throughout the discussion, we have assumed that the 
distance to Tycho's SNR is 2.1 kpc. The measurement 
is based on the estimated proton temperature (from the 
observed line width of the broad component) and the 
optical proper motion, assuming no cosmic ray accelera- 
tion at the shock. The cosmic ray acceleration, if signif- 
ican t, can effectively r educe the postshock temperature 
(see IHelder et al.|[2009l and references therein); thus the 
distance may have been underestimated. However, as 
has been discussed above, the optical observations are 
consistent with cosmic ray acceleration not being very 
efficient in this region, so the distance of 2.1 kpc, derived 
assuming no cosmic ray acceleration, remains a reason- 
able value. 

Using Chandra observations, iWarren et al.l ()2005[ ) 
found the locations of the shock front (SF) and the con- 
tact discontinuity (CD) along the boundary of the rem- 
nant. They interpreted the small separation between the 
two as an indication of efficient cosmic ray acceleration. 
The region of knot g is one of the regions where the SF- 
CD separation is smallest (except those regions of ejecta 
protrusion). A simple extrapolation of their argument 
will lead to the most efficient cosmic ray acceleration 
in this region, being inconsistent with our results. How- 
ever, the region around knot g is where the remnant could 
be iii t eracting with d ense ambient clouds (iRevnoso et al.l 
ll997HLee etalll200 4). Thus, the small SF-CD separa- 
tion in this region could be due to a recent encounter of 
the shock with the dense ambient gas, instead of efficient 
cosmic ray acceleration. 

In conclusion, we have presented high resolution Ha 
imaging observations of Tycho, revealing the existence of 
a thin precursor which we interpret as a cosmic ray pre- 
cursor. While the current observation is consistent with 
inefficient acceleration, the observation of the precursor 
itself provides an important opportunity to constrain the 
key parameters of the acceleration, such as the diffusion 
coefficient and the injection parameters. A comparison 
with detailed numerical simulations will be critical to 
study the detailed physics of cosmic ray acceleration. 

This research was supported by STScI grant GO- 
11184.01-A-R to the Smithsonian Astrophysical Obser- 
vatory. 
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